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Abstract 

The usefulness of the three-body mechanism of ry-meson produc- 
tion is tested by extending to high energies a recently developed two 
step model of the pd — > 3 He n reaction by us , which was successful in 
reproducing the energy dependence of the data near threshold. The n 
production in this model, proceeds in two steps via the pp — ► ird and 
ttN — * rjp reactions with the intermediate particles being off-shell. 
The final state interaction(FSI) is incorporated through a T-matrix 
for n 3 He elastic scattering constructed using few body equations. In 
a comparison of the calculated angular distributions with some re- 
cently available data, we discover the limitations of the model at high 
energies. 

PACS numbers: 13.75.-n, 25.40.Ve, 25.10,+s 
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1 Introduction 



The importance of the three-body mechanism (3-bm) in meson production 
was first shown in [1, 2] by Laget and Lecolley. Whereas in the pion pro- 
ducing reactions pd -h> 3 H n + and pd — > 3 He n°, this mechanism accounted 
for the discrepancies left out at backward angles by the one- and two-body 
mechanisms, in the case of the pd — > 3 He rj reaction, it was found to play 
the dominant role in reproducing the Saturne data at threshold [3] and that 
at high energies [4] for 6 V = 180°. The one- and two-body mechanisms in 
the pd ^ 3 He r\ reaction, were found to underestimate the experimental cross 
sections by two orders of magnitude. Motivated by these works, more re- 
cently, two step models [5, 6, 7] which essentially involved the 3-bm, were 
constructed to study the pd — > 3 He rj data near threshold [3, 8]. The r^-meson 
in these models is produced in two steps via the pp — > ivd and nN — > r/N re- 
actions, thus involving the participation and sharing of momentum by three 
nucleons which is typical of a 3-bm (see Fig. 1). The calculation in [7] is the 
most involved one available on the pd — > 3 He rj reaction in literature. In it, (i) 
the final state 7? 3 He interaction which is crucial for the energy dependence of 
this reaction is incorporated through a T-matrix constructed using few-body 
equations and (ii) the fermi motion and off-shell nature of the intermediate 
particles in the 3-bm is taken into account properly. Some other theoretical 
investigations of the pd — > 3 He i] reaction at threshold and high energies can 
also be found in [9, 10]. The work in [10] needs a special comment to which 
we shall return later. 

The explanation for the dominance of the 3-bm in rj production is 
straightforward. Due to the large mass of the r\- meson (~ 547 MeV), the 
momentum transfer involved in this reaction is as large as 900 MeV/c near 
threshold. The 3-bm allows this large momentum to be shared among three 
nucleons, thus making it easily digestible by the nucleus than in the one- 
and two-body mechanisms where the nuclear form factor (evaluated at large 
momenta) is rather small. Since the two step model alone could reproduce 
the strong energy dependence of the threshold data [7] and given the fact 
that the momentum transfer in the pd — > 3 He r] reaction continues to be large 
(~ 650 MeV/c) even at 1 GeV above threshold, one would speculate that 
the same model would also be successful in reproducing the recently reported 
angular distributions [11] at high energies. An extension of the model in [7] 
to compare with the high energy data of [11], therefore, seems timely. 



2 




\.p p -> K d 
2. ti n -> T) p 



J He 



1 . p p — > 71 a 

2. 71 p -> f\ p- 



He 



Figure 1: The diagram corresponding to the three-body mechanism of rj- 
meson production in the pd — » 3 He rj reaction. 



However, before proceeding further, let us recapitulate the following 
observations from the different 3-bm calculations carried out in the literature 
at energies beyond threshold: (a) the success of the 3-bm reported in [2] for 
the the pd — > 3 He rj reaction over a range of proton beam energies up to 2.5 
GeV, was confined to one particular angle, namely, 9 V = 180°, (b) in case of 
pion production at high energies, the 3-bm was mostly responsible in remov- 
ing the discrepancies between theory [1] and data at backward pion angles, 
(c) the theoretical estimates [6, 12] of angular distributions (at threshold and 
at high energies) of the pd K + reactions made within a two step model, 
seem to indicate backward peaked cross sections and (d) the cross sections 
calculated in [7] for the pd — > 3 He rj reaction within the two step model were 
nearly isotropic near threshold; however, there was a hint of a backward peak 
as one moved about 10 MeV away from threshold. 

That is, in all, the existing 3-bm calculations point toward a back- 
ward angle dominance of the cross sections in this model. In contrast to this, 
the measured angular distributions of the pd — > 3 He rj reaction from about 
90-200 MeV above threshold are all forward peaked [11, 13]. Thus, a priori, 
one would doubt the possibility of the two-step model being able to repro- 
duce the high energy data. Indeed in [14], in an attempt to reproduce the 
differential cross section data at beam energies of 930, 965, 1037 and 1100 
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MeV, using a Monte-Carlo description of the pd — > 3 He rj reaction within a 
two-step model, the authors find backward peaked cross sections which are 
in complete disagreement with data, but are however, in line with the above 
theoretical findings. A very recent work [10], which reports the success of the 
two step model in reproducing the data on angular distributions mentioned 
above, is then isolated in its findings. The results in [10] have been obtained 
using a factorization approximation, which neglects the off-shell nature of 
the intermediate pion propagator as well as the t-matrices for the %N — > rjN 
and pp — > n + d reactions. It restricts the pion in an ad hoc way to go at 
0°. Besides, though not mentioned clearly, it seems that the final state in- 
teraction has also been incorporated within a factorization approximation as 
in [5]. The model with such approximations, does not truly represent the 
three-body Feynman diagram of Fig. 1 (where the intermediate particles 
are off-shell), hence the results of such specious models should be treated 
with caution. To obtain a proper estimate of the contribution of the 3-bm 
at higher energies, we extend our earlier calculation [7] to evaluate the dif- 
ferential cross sections at 930, 965, 980, 1037 and 1100 MeV beam energy 
for which data are available. These calculations include the off-shell and fi- 
nal state interaction effects elaborately and do not restrict the pion to any 
specific angle. Unlike [10] and as expected, partly from our earlier work [7] 
and partly from calculations of other reactions with 3-bm in literature (men- 
tioned in (c) above), we get backward peaked angular distributions which are 
at variance with the corresponding forward peaked data. The cause for this 
discrepancy may lie in the interaction among the particles in the intermedi- 
ate state or the increasing importance of some other mechanisms which need 
to be investigated. We therefore conclude that the 3-bm calculations of the 
pd ^ 3 He rj reaction at energies away from the threshold do not reproduce the 
observed forward peaking of the rj angular distribution seen in the Uppsala 
data. This is in contradiction with the observations reported in [10]. 

In the next section we describe the inputs of the model only briefly 
and move on to the discussion of the results. The details of the model can 
be found in [7]. 
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2 Two step model of the pd -* 3 Her] reaction 

In the two step model of rj production, the incident proton interacts with 
one of the nucleons in the deuteron (d) to produce a deuteron {d') and pion, 
which in turn interacts with the other nucleon in d to produce an ?7-meson 
and a proton. This proton and the d' combine to form the 3 He nucleus. The 
transition matrix for the pd — > 3 He r\ reaction which includes the final state 
r] 3 B.e interaction is given as, 

T = < k v ;m 3 \ T pd ^3 Her) \k p ;m 1 m 2 > + (1) 

r dq < k v ; m 3 \T v 3 He \ q; m' 3 > , -* 

E J (27r) 3 E{k ) _ E{q) + ie <r,m 3 \T pdMv | k p ; mi m 2 >. 

m 3 

The matrix elements < \T pd _^3 Her) \ > in the above equation correspond to the 
Born amplitude for the pd 3 He r\ reaction. k p and k v are the asymptotic 
momenta of the particles in the initial and final states and mi, m 2 and m 3 
are the spin projections of the proton, deuteron and 3 He respectively. The 
T-matrix for ?7 3 He elastic scattering, T v 3 He , is evaluated using four-body 
equations for the r/(3N) system. The Born amplitude is given within the two 
step model as, 

< \T pd _3 Herj \ >= i J ^ ^ Yl <pn\d><nd\T pp ^ wd \pp> (2) 

V / V / intra' s 

1 „ 

x 7l2 2 ■ - \ K VpITttN^pW N > < ' He\pd> , 

[K w — m n + ie) 

where the sum runs over the spin projections of the intermediate off-shell 
particles. The four momentum of the intermediate pion, k n , which could 
either be a 7r + or 7r° is given as, 

k 2 ^ 

K = y + -k^+P! +p 2 (3) 

where p\ and p 2 are the fermi momenta of the intermediate nucleons. Each 
of the individual matrix elements in the above equation is expressed in terms 
of partial wave expansions. The matrix elements for the pp — > nd reaction, 
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parametrized in terms of the available experimental data are taken from 
Ref. [15]. For the nN — > rjp reaction, we use the single resonance, coupled 
channel model of Ref. [16]. This model treats the nN, rjN and tttcN as 
coupled channels. We use the T-matrix in this model with the inclusion of 
the Su, Pu and TJ 13 resonances. Since the r] — Pu coupling used in Ref. [16], 
as noticed in [10], does not reproduce the data on n N — > rj N reaction, we 
use it as given in [10]. The matrix elements < pn\d > and < 3 He|jo<i > consist 
of the deuteron and helium wave functions in momentum space which depend 
on the momenta p\ and P2 respectively. The details of the wave functions and 
the r] 3 He T-matrix are given in [7]. In the calculation of the ?7 3 He T-matrix, 
we use the coupled channel elementary t-matrix, t^N^N as in [17]. This 
t-matrix was the one which gave the best agreement in [7] with the threshold 
data. 



3 Results and Discussion 

Recent measurements of the p d — > 3 He rj reaction, from about 22 MeV to 1 15 
MeV above threshold [11, 13] show that, as we go away from the threshold, 
the angular distribution starts developing anisotropy, which increases with 
energy. The GEM data [13] show a strong forward peak and the data from the 
Wasa-Promice collaboration [11] are found to be maximal for cos 6 V ^+0.5. 
At threshold, the angular distribution is isotropic. 

Fig. 2 shows our calculated differential cross sections, using the two 
step model described in the previous section, along with the experimental 
data from [11, 13]. The dashed curves represent the plane wave calculations, 
while the solid lines correspond to the calculations including the r; 3 He final 
state interaction (FSI). We find that the observed trend in the anisotropy of 
the data is not at all reflected in the calculated results. Unlike the data, the 
latter peak at backward angles. This in contrast to the findings in [10]. 

Fig. 3 shows the data and calculation of the spin averaged amplitude 
which is defined as, 

I n |2 _ K °tot / t \ 

1/1 k v 47r 1 ' 

where k p and k v are the proton and rj momentum in the center of mass 
system. a tot is the angle integrated cross section. In the experiment it is 
obtained by making a polynomial fit to the data at the available angles 
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Figure 2: The angular distribution of the pd — > 3 He?7 reaction. The data 
are from Ref. [11, 13]. The data points (A) and (B) of [11], represent respec- 
tively, the data taken with only the 3 He detected and with 3 He detected in 
coincidence with two photons. The dashed lines (solid lines) are calculations 
without (with) inclusion of the final state interaction (FSI). 



and then integrating the parameterized distribution. We observe two things: 

(i) the two step model results are in excellent accord with the data near 
threshold. The model however underestimates the data at higher energies. 

(ii) Comparison of the calculated results with and without FSI show that, 
while near threshold the FSI increases the cross section, at higher energies 
its effect is to decrease the cross sections. This, as discussed in our earlier 
work [7] , is the reflection of the decreasing importance of the off-shell rj— 3 He 
scattering in the final state at high energies. 

To see the implication of the restriction on the direction of the inter- 
mediate pion, we show in Fig. 4, our calculated results at one beam energy, 
for pions restricted up to various emission angles. As found in [10], the 
forward pion emission indeed skews the differential cross section to forward 
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Figure 3: The square of pd — > 3 Her? amplitude as defined in eq. (4). The 
dashed (solid) lines correspond to the calculations without (with) FSI. The 
data are from refs [8, 11, 13]. 




Figure 4: The angular distribution for the pd — > 3 He rj reaction at 930 MeV 
beam energy. The solid line shows the cross section with the intermediate 
pion angle restricted from to 10° with its scale on the right hand side. The 
long dashed, dot-dashed and dotted curves show the cross sections obtained 
after restricting the pion angle from to 30°, to 50° and allowing the full 
range of to 180° respectively (with the scale given on the left). 
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angles. As such, this restriction is of course ad hoc. Any legitimacy in this 
may be sought in the interaction among the particles in the intermediate 
state. We also see that the cross sections fall with reduction in the range of 
allowed pion angles. The scale given on the right side of Fig. 4 for the to 10° 
restriction (solid line) shows that though the shape of the experimental an- 
gular distribution is correctly reproduced, the magnitude is underestimated 
by about a factor of 40. This is again in contrast to the finding in [10], where 
the agreement with absolute cross sections, with the pion angle restricted to 
0° is not so bad. The reason for this discrepancy could lie in the neglect of 
the off-shell effects at the intermediate vertices in the reaction mechanism in 
[10]. Since we only wish to demonstrate in Fig. 4 the effect of restricting 
pion angles, these calculations have been done without including the FSI. 

4 Conclusions 

The role of the three-body mechanism had been demonstrated unambigu- 
ously in our earlier work on the pd — * 3 He rj reaction near threshold [7]. 
Indeed, isotropic angular distributions and the magnitude of the total cross 
sections were well reproduced. Extending the model to higher energies (with- 
out making any rough approximations outside the spirit of the model) leads 
to a large discrepancy between model predictions and available data. We 
believe that any ad hoc manipulations of the model [10], like e.g., restricting 
the pion to be on-shell and to be produced in the forward direction (9 pn = 0), 
do not represent the three-body mechanism in the context of the Feynman 
diagram and hence are hard to justify. Had these conditions really been the 
preferred ones in the three-body mechanism, it should have come out anyway 
through a proper calculation as performed in the present work. However, this 
is not the case. 

The failure of the model, especially in reproducing data at forward 
angles hints toward some missing component in the model. In the present 
work, we do not take into account the interaction of the intermediate off- 
shell particles which in principle could be quite important. Contribution of 
coupled channel processes such as the pd ^ 3 He 7r° may also turn out to 
be significant. The possibility of other mechanisms, such as the one- and 
two-body mechanisms taking over at high energies seems unlikely for the 
following reasons: (i) the momentum transfer in the pd — > 3 He r\ reaction, 
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as mentioned in the introduction, continues to be quite large even at high 
energies (ii) as far as kinematics is concerned, the pd — > 3 He rj and pd — >aH 
K + reactions are similar. Hence, the results of the theoretical investigation 
of the pd ^jji K + reaction up to proton beam energy of 3 GeV [12], where 
the authors find the one- and two-body mechanisms to contribute 2-3 orders 
of magnitude lesser than the three-body mechanism at all angles, taken along 
with similar findings of [2] f or pd — > 3 He r/ at 9 V = 180°, any expectations from 
the one- and two-body diagrams for r\ production in reproducing the high 
energy angular distribution data for the pd — > 3 He rj reaction seem remote. 
Though the effect of the final state interaction has not been studied with 
the one- and two-body mechanisms, from the present work it seems that at 
high energies this will only further reduce the cross sections (see Fig. 3 for 
example), making the one and two-body mechanisms even more negligible. 

The above conclusions could probably change if the deuteron and 
the 3 He have some hitherto unconsidered dynamics (like quarks) at short 
distances, which influence their structures. Our efforts in future would be 
to explore such avenues and the role of other possible diagrams such as the 
direct production of the r\ meson. 
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